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POLYHYPERBOLIC CARDINAL SPLINES 


JEFF LEDFORD 

Abstract. In this note we discuss solutions of differential equa¬ 
tion — a^)^u = 0 on M \ Z, which we call hyperbolic splines. 
We develop the fundamental function of interpolation and prove 
various properties related to these splines. 

Keywords, cardinal interpolation, hyperbolic splines 


1. Introduction 

Hyperbolic splines have been studied extensively, and are often called 
exponential splines [1]. In this article, we prove results similar to those 
found in n. 0. and [3] . In each of these papers, a fundamental function 
of interpolation is studied. 

Our main result concerns band-limited functions and shows that as 
the degree of tends to inhnity the polyhyperbolic spline interpolant 
tends to the target function. Our methods require the use of the Fourier 
transform, for which we set the following convention. If g & Li(R), then 
its Fourier transform, written g, is given by 

m = f s(x)e-'‘^d(. 

Jr 

This dehnition is extended to distributions when necessary. 

2. Polyhyperbolic Cardinal Splines 

Definition 1. Given a > 0 and A; G N, the set of all k-hyperbolic 
cardinal splines is given by 

H^= [f e : {D^ - = 0, on R \ Z} (1) 

If A: > 1, we say that / G is a polyhyperbolic cardinal spline. The 
justihcation for the name comes from the fact that {cosh(a-), sinh(a-)} 
forms a fundamental set of solutions for the differential equation {D^ — 
= 0. In fact, it is not too hard to see that : 0 < j < A;} 

forms a fundamental set of solutions for the equation {D^ — a‘^)^f = 0. 

Most of this paper is analogous to the treatment given to polyhar¬ 
monic cardinal splines in [2]. We begin with what is called the funda¬ 
mental solution in [2]. The solution to the equation — a‘^)^Ek{t) = 
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6(t) is given by 

Ek{t) = ( 2 ) 

where f*^ is the convolution of / with itself k times and Ca,k is chosen 
so that the resulting multiple of the Dirac h-function is 1. The Fourier 
transform of this function plays an important role in what follows, it 
is given by 

= + ( 3 ) 

From this function we develop a fundamental function of cardinal 
interpolation L^, which is dehned by its Fourier transform 


LK) 


(2ir) 


1/2 __ 

Xl/ez Ek{^ — Svrj) 


( 4 ) 


We will also be concerned with the related periodic distribution <hfc, 
dehned by 

MO = mWO- ( 5 ) 

In order to study Lk and <Ffc further, we make use of some basic analysis 
in the complex domain. For A C M, we dehne 


A, = {C e C : 3fJ(C) e A and |^(C)| < e}. 

We may now prove our hrst lemma. 

Lemma 1. There exists e > 0 such that both CLi^d Lk have extensions 
to Me which are analytic. 


Proof. We follow along the same lines as the a = 0 case which may 
be found in |2]. We will henceforth assume that a > 0. Set q{0 = 
— + a^), where C = ^ + O- We have 

(27r)-'/2|5;Kr' = {1 + i9K)i‘D0} 

where 

F(C) = ^l?(C - 2irj)]-‘. 

jVo 

Since a > 0, q~^ is analytic in Me for any e G (0, a). It is also clear that 
F is analytic on (—vr, 7r]e. Now since > 0 for G (—vr, vr], we 

can hnd e > 0 (which may be smaller than our hrst choice) such that 
1 + [g((C)]^F(C) has no zeros in (—7r,7r]e. This shows that <I>fc possesses 
an extension which is analytic in (—7r,7r]e. Using periodicity, we may 
extend this result to Me. The analytic extension of Lk to Me is given 
by Lk{0 = 4>fc(C)[5'(C)]~^) which as the product of analytic functions 
is analytic. □ 
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Lemma 2. We have that 

^k{x) = ^aj5{x - i), ( 6 ) 

where the aj’s depend on k and a. Furthermore, there exists two con¬ 
stants c,C>0 such that for all j G Z, 

kil < (7) 

Proof. Periodicity implies that ([6]) holds with 

a, = (27r)-'/2 f 

J (—7r,7r] 

To see that ([7]) holds, we use Lemma 1 to replace (—tt, tt] with G 
C : 3?(C) G (—TT, tt] and 7y(C) = sgn(j)c}, where 0 < c < e. We have 

\aj\ < (27r)-i/2g-c|i| f 

J{-n,TT] 

□ 

Proposition 1. Let be defined by its Fourier transform (jl]), where 
A; G N. Then has the following properties: 

(i) Lfc is a k-hyperbolic cardinal spline. 

(ii) For all j G Z, Lk[j) = doj, where 6 is the Kronecker 6 function. 

(iii) There are constants c,C>0 (which depend on k and a) such 
that 

\Lk{x)\<Ce-<^\^\ 

for all a; G M. 

(iv) Lfc has the following representation in terms of E^: 

Lk{x) = ^ ajE{x - j) = * Ek{x), 

j&i. 

where <hA; is the function whose Fourier transform is defined by 
and the Oj’s are those from Lemma 2. 

Proof, (i) We see that {D'^ — a^)^Lk = from examining the equation 
in the Fourier domain. Now Lemma 2 shows that L^ G H^. 

(ii) We have 
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The interchange of the snm and in the integral in the “periodization 
trick” is justihed by the decay of Ek and the periodicity of the denom¬ 
inator of Lfc. That the sum is (27r)“^/^ is clear from (|1]) and Lemma 
1 . 

(hi) An argument completely analogous to the one given in Lemma 
2 show that this is true. 

(iv) This is clear from (|1]) and ([S]). □ 

We turn briefly to the data that we wish to interpolate. 

Definition 2. For > 0, let be the set of all sequences b = {bj : 
j G Z} which satisfy 

|4yl<C(l + |i|)'’, (8) 

for some constant C > 0. 

Our next result shows what to expect when we interpolate data of 
this type. 

Proposition 2. Suppose that b = {bj : j ^ 7 j} is a sequence of polyno¬ 
mial growth and define the function f^ by 

fb{x) = ^bjLk{x - j). (9) 

The following are true: 

(i) The expansion (jl]) converges absolutely and uniformly in every 
compact subset o/M. 

(ii) The function fi is a k-hyperholic spline and fb{j) = bj for all 
J e Z. 

(hi) If b E Yh, then fi satisfies 

\Mx)\<C(l + \x\f, 

for some C > 0. 

(iv) If g E span({(-)'^e^“(') : 0 < j < k}) and b = {g{j) : j E Z}, 
then fb{x) = g{x) for all x G M. 

Proof, (i) and (ii) are straightforward consequences of Proposition 1. 

(iii) Let E Z satisfy < x < + 1. We have 

lAWI < E - j)\ < 

i + b1 


= 0(1 + 




1 -|- IrUx 


^c\x-m^\-c\x-j\ 


=:0(1 + |m,|)V"l™l{l + A + B}, 
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□ 


where 


and 



/3 

^c\x-ma;\-c\x-j\ 



^c\x-m:^\-c\x-j\ 


The estimates being similar, we show only the estimate for A: 




OO 

E 

J='ma;+l 


1 + b'l 

1 + Im^ 


^-c\x-j\ 



1 + \mx + k\ 
1 + \m^\ 


^-c(k+mx-x) 


k=l 


<e2"^(l + |A;|f < oo. 

k=l 


Here we have used the triangle inequality and reindexed the sum. In the 
second inequality we have used that /3 > 0. Finally, the last inequality 
comes from the fact that rux < x < rux + 1. Combining this estimate 
with the analogous estimate for B yields the desired result. 

(iv) This is a consequence of the linear independence of the set 
: 0 < j < /c} C H^. Since fb{j) and g{j) conincide on 
2k consecutive integers, say 0 < j < 2A; — 1, we have fb{x) = g{x) 
on the interval [0, 2A; — 1]. Repeating this argument in both directions 
shows that fb{x) = g{x) for all a; G M. 

We end this note with a result concerning the behavior of polyhy- 
perbolic splines as their degree tends to inhnity. 

For this we recall that the classical Paley-Wiener space PW, is de- 
hned by 


PW = {ge L 2 (M) : ^(0 = 0 a.e. ^ ^ [-7r,7r]}. 

A function g G PW is often called hand-limited. The Paley-Wiener 
theorem shows that band-limited functions are entire functions of ex¬ 
ponential type whose restrictions to the real line are in L 2 {M.). 
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In light of Proposition 2, we may interpolate a function g G PW 
with a fc-hyperbolic spline as follows 

h[g]{.x) = ^g{])Lk{x - i). (10) 


Theorem 1. Suppose that g G PW. The function Ik[g\ defined in ffTOj) 
satisfies 

(i) 

IIS'- 4[5']IU2(R) = 0, and (11) 


lim \g{x) — Ik[g]{x)\ = 0, uniformly on M. (12) 

/c—>-oo 

Proof. We begin by noting that if y = {pj : j G Z} G for 1 < p < cxd 
then fy G Lp(R), where fy is dehned by ([9]). Since Lk has exponential 
decay, the estimates of the form 

||/y||Li(R) <Ck\\y\U„ and 
||/t/IUoo(R) <Cl\\y\\i^ 

are easily established. Now the Riesz-Thorin interpolation theorem 
provides a constant Ck,p > 0 such that 

||/y||Lp(R) < Ck,p\\y\\ep- 

It is well known that if p G PW, then {g{j) : j G Z} G £ 2 - Hence we 
can see that 


||4[5']IU2(R) < Ck\\{g{j) : j G Zjll^j 

for some constant Ck > 0. This allows us to make use of the Fourier 
transform. We have 


4[9lK) = i*K)5^9(j)e«. 

jez 

In particular, for ^ G (—vr, tt], Ik[g]{^) = V^Lk{^)g{^). Now we have 
the following estimates using Plancherel’s identity and the fact that g 
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is band-limited: 



The second eqnality follows from periodization, while the third follows 
from Tonelli’s theorem and manipnlating (jlj). The hnal ineqnality fol¬ 
lows from the positivity of the terms in the snm. In light of (l3i). we 
have the following elementary estimate for ^ G [—vr, vr] and i ^ 0: 


Lk{^-27ri)\<{27r)-^/^ 




{2\i\ - l^TT^ + a^ 


k 


This allows ns to nse the dominated convergence theorem since the 
corresponding series is bonnded. For |,^| < tt, the estimate above works 
with replacing vr^ in the nnmerator, thns each term in parentheses 
is strictly less that 1, hence tends to 0 as fc —)■ cxd. This completes 
the proof of (i). The proof of (ii) follows from (i) together with the 
Canchy-Schwarz ineqnality. Using the inversion formnla and the trian¬ 
gle ineqnality, we have 


\g{x) - Ik[g]{x)\ = (27r) 


- 1/2 


m - 


<(27r)-V2/ ^(0-4[^](0 


dC 


Now we may nse periodization and the Canchy-Schwarz ineqnality and 
reasoning similar to that nsed in (i) to see that 


\g{x) - Ik[g]{x)\ <2 I ^ - 27r4 j d^, 

.£7^0 


which tends to 0 independently of x G M. 


□ 
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